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Abstract 

We study multiple orthogonal polynomials of type I and type II which have ortho- 
gonality conditions with respect to r measures. These polynomials are connected by 
their recurrence relation of order r + 1 . First we show a relation with the eigenvalue 
problem of a banded lower Hessenberg matrix L n , containing the recurrence coef- 
ficients. As a consequence, we easily find that the multiple orthogonal polynomials 
of type I and type II satisfy a generalized Christoffel-Darboux identity. Furthermore, 
we explain the notion of multiple Gaussian quadrature (for proper multi-indices), 
which is an extension of the theory of Gaussian quadrature for orthogonal polyno- 
mials and was introduced by C. F. Borges. In particular we show that the quadrature 
points and quadrature weights can be expressed in terms of the eigenvalue problem 
of L n . 
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1 Introduction 



Multiple orthogonal polynomials arise naturally in the theory of simultaneous rational 
approximation, in particular in the Hermite-Pade approximation of a system of r £ N 
(Markov and Stieltjes) functions [7,8,14,15]. They are a generalization of orthogonal poly- 
nomials in the sense that they satisfy orthogonality conditions with respect to r £ N 
measures fii, . . . , \i r for which all the moments exist. In the literature one can find already 
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a lot of examples of such polynomials [1-4,17]. Normally the measures are taken to be 
positive. However, in this paper we study formal multiple orthogonal polynomials which 
means that we allow complex measures. 

We will only consider multiple orthogonal polynomials with respect to proper multi- 
indices. The proper multi-index corresponding to n G N = N U {0} is 

V n — (to + 1, m + 1, . . . , m + 1, to, to, . . . , to) G Nq, 



s times 



r— s times 



where n = mr + s, < s < r. There exist two types of multiple orthogonal polynomials, 
type I and type II. Let F±, . . . , F r be the supports of the r measures. A multiple orthogonal 
polynomial P n = Pp n of type II with respect to the proper multi-index v n , is a polynomial 
of degree at most n, which satisfies the orthogonality conditions 



/ P n (x) x E dnj(x) = 0, £ = 0,...,u n (J)-l, 



j = l,...,r. 



Here v n {j) is the jth component of V n . Equation (1.1) leads to a system of n homogeneous 
linear equations for the n + 1 unknown coefficients of P n . A basic requirement to have a 
good definition is that every possible solution of the system (1.1) has exactly degree n. This 
is equal to the assumption that (1.1) has a unique solution (up to a scalar multiplicative 
constant) which has exactly degree n. In that case we call V n a normal multi-index for 
III, jJL r . Let rrSp = J r . x e dfij(x) be the £th moment of the measure fij and set 



D r . 



where 



D 



(i) 



D 



(r) 

n,i>n(r) 



(1.2) 



D 



U) 0) 



m £-l 



m 



is an n x £ matrix of moments of the measure Then D n is the matrix of the linear 
system (1.1), without the last column. It is known and easily verified that V n is normal if 
and only if D n has rank n [3,4,15]. In the case that all the proper multi-indices v ni n G N, 
are normal, we call the system of measures a weakly complete system. 



,A nr ) T , corres- 



A type I multiple orthogonal vector polynomial A n = A$ n = (A njl , 
ponding to the proper multi-index u n , consists of r polynomials A n j of degree at most 
v n (j) — 1, satisfying the orthogonality conditions 



x £ J2 A n,j(x) dfi j (x) = 0, £ = 0, 1, . . . , n - 2. 

J .7 = 1 



(1.3) 
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The usual requirement is that every A n j has exactly degree v n (j) — 1. However, as men- 
tioned before, in this paper we only consider proper multi-indices. So, if the system (1.3) 
has a unique solution up to a multiplicative scalar constant and A n ^ s has exactly degree 
m, where n = mr + s, < s < r, we obtain a good definition. Note that these conditions 
are satisfied if z? n -i is normal, which means that D n _i has rank n — 1. This easily follows 
from the fact that the matrix of the linear system (1.3), without the s(m + l)th column, 
is equal to D^_ v 

In Section 3 we approximate integrals of the kind 



simultaneously, using weighted quadrature formulas which have the same n quadrature 
nodes. C. F. Borges studied this already in [5]. The point of interest is then to find the 
appropriate quadrature nodes and weights maximizing the vector order of the weighted 
quadrature formulas along the proper multi-indices. In [5, Theorem 2], C. F. Borges 
showed that this is obtained if we require each weighted quadrature formula to be in- 
terpolating and the quadrature nodes to be the zeros of the type II multiple orthogonal 
polynomial P n , assumed to be simple. We recall this notion of multiple Gaussian quadra- 
ture in Theorem 3.1. 

In the case of Gaussian quadrature (r = 1) the nodes and weights can be expressed in 
terms of the eigenvalue problem of a Jacobi matrix, containing the recurrence coefficients, 
see e.g. [10, Chapter 3, §2.3] and [11]. In Theorem 3.2 we extend this to multiple Gaussian 
quadrature. Here we need some properties of multiple orthogonal polynomials. In Section 
2 we recall that, for a weakly complete system, the multiple orthogonal polynomials of 
type I and type II each satisfy a recurrence relation of order r + 1 which are closely 
connected, see, e.g., [14, §24] and [6,12,16,18]. As a consequence, these polynomials of 
type I and II are then linked to the left and right eigenvalue problem, respectively, of 
a banded lower Hessenberg matrix L n containing the recurrence coefficients as in (2.8). 
This then also leads to a generalized Christoffel-Darboux identity in Theorem 2.5, similar 
to [16, (21)]. Using all this, we show that in the case of multiple Gaussian quadrature 
the nodes are the eigenvalues of L n and the weights can be expressed in terms of the 
corresponding left and right eigenvectors. 



2 Recurrence relation of order r + 1 

2. 1 Relation between type I and type II 

Suppose that the system of measures fj,i, . . . , /i r forms a weakly complete system. The 
monic multiple orthogonal polynomials P n of type II corresponding to proper multi- 
indices are then uniquely determined. Furthermore, as mentioned in the introduction, 
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also the multiple orthogonal polynomials of type I, A n , are unique up to a normalizing 
multiplicative constant. In this paper we take the normalization 

V 

fx n - ly £Anj(x)d f JL j (x) = l. (2.1) 
J 3=1 

Note that the integral in (2.1) cannot be zero since z7 n is a normal multi-index. 

It is known that multiple orthogonal polynomials satisfy a recurrence relation of order 
r + 1, see, e.g., [14, §24] and [6,12,16,18]. For the monic multiple orthogonal polynomials 
of type II this is 

r 

xP n (x) = P n+1 (x) + anj P n -j{x), n > 0, (2.2) 

3=0 

with initial conditions Pq = 1 and Pj = 0, j — — 1, —2, . . . , —v. The proof is similar as in 
the case of orthogonal polynomials. Furthermore, there exist some integral representations 
for the recurrence coefficients in terms of the multiple orthogonal polynomials of type I 
and type II. If we integrate (2.2) with respect to the measures Y?i=\ A n+ i_ k ,e( x ) d/i^x), 
k — 0, 1, . . . , r, then we find 

/ xP n (x)Y, A n+i-kA x ) d^e(x) = / P n+1 (x) A n+i-k,e( x ) dfie(x) 
i=i J i=\ 

r - r 



j=0 1=1 



Now apply the orthogonality relations of the multiple orthogonal polynomials of type I 
and II. Because of the normalization (2.1) and the fact that the P n are monic, we then 
obtain 

a nd = / xP n (x) Y A n+ i_ j>e (x) dfjL t (x), j = 0, 1, . . . , min(r, n). (2.3) 
J i=\ 



The multiple orthogonal polynomials of type I also satisfy a recurrence relation of order 
r + 1 which is related to the one of type II. Denote by C r [x] the space of r-dimensional 

— * 

vector polynomials with complex coefficients. The vector polynomials A n then form a 
basis for C r [x], so that there exist unique c n j for which xA n (x) = J2k=i °n,k Ak{x). By 
linearity we then have 

/r n+r „ r 

xPi(x) Y A n ,e( x ) dni(x) = Y °n,k / Pi(x) Y A k /{x) dfi e {x), i = 0, 1, . . . 
i=\ k=i J e=i 

The cases % = 0, 1, . . . , n — 3 give rise to c n> i = ■ ■ ■ = c ni „„2 = by the orthogonality 
relations (1.3). Furthermore, use the orthogonality relations (1.1) and the normalization 
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(2.1) to get 

C n , ra _i = 1, C n>n+j = / xPn+j-^^An^x) dfl e {x) , j = 0, 1, . . . , T. 

J 1=1 

Comparing this with (2.3) we obtain c n>n+ j = a n+ j-ij, j = 0,1, ... ,r. So if the monic 
multiple orthogonal polynomials of type II satisfy the recurrence relation (2.2) then the 
multiple orthogonal polynomials of type I with normalization (2.1) satisfy the recurrence 
relation 

r 

xA n (x) = A n ^i(x) + On+j-lj A n+j(x), 71 > 1, (2.4) 

3=0 

with initial conditions Aq = and A\,...,A r . Denote by the matrix obtained 

from D n by deleting the jth row and the zth column (and set det-Df 1 ' 1 ^ equal to 1). For 
1 < i, j < r we then have 



Aij — < 



det D (j,i) 

[ ' det A ' 3 ~ h (2.5) 
, j > i, 



which are functions of the (first r) moments of the measures fj,i, . . . , /i r . Notice that 
A?'J 7^ 0j i = 1) • • • > r > because the multi-indices V n are normal. 

Remark 2.1 Let n = mr + s, < s < r. Applying the orthogonality conditions of the 
polynomial P n+r _i, we obtain from (2.3) 



'•n+r- 



-i,r = J xP n+r ^i(x)Y^ A n/ (x) dfi e (x) = J xP n+r _i(x)A ntS (x) dfi s (x), n>l. 



We assumed that V n _ x and normal, so A nyS has exactly degree m and 

a n+r _i ir ^ 0, n > 1. (2.6) 

This means that we can recover the vector-polynomials A n by (2.4) if we know the initial 
conditions (2.5) and the recurrence coefficients. 

Remark 2.2 We supposed that the measures . . . , /j, r form a weakly complete system. 
If we replace these measures by 

r 

«1,1 /"l, «2,1 + «2,2 1^2, • • • , a r,3 Vji ( 2 -7) 

where a it j G C, 1 < j < i < r, then it is clear that this system of measures forms a weakly 
complete system if and only if ajj ^ 0, j — 1, . . . , r. Furthermore, all these systems of 
measures have the same multiple orthogonal polynomials of type II with proper multi- 
indices and so the same recurrence coefficients. It is obvious that every set of measures of 
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the form (2.7) then corresponds to another choice of initial conditions for the recurrence 
relation (2.4) where Ajj ^ 0, j — 1, . . . , r. 



2.2 Eigenvalue problem of the banded Hessenberg matrix L n 



In this section we introduce the banded lower Hessenberg matrix 



ao,o 1 

ai,i ai )0 1 

a 2,2 a 2,l ^2,0 1 

a r r a r i a r o 1 

a r+ i jr ttr+1,1 ^r+1,0 1 











1 



\ 



V 



o 



d n -l,r On-1,1 On-1,0 j 



(2.8) 



which contains the recurrence coefficients of the multiple orthogonal polynomials. If we 
expand the determinant det(xl n — L n ) along the last row, then we find that these polyno- 
mials satisfy the recurrence relation (2.2). So P n (x) = det (xl n — L n ) and the eigenvalues 
. . . ,x n>n of the matrix L n then coincide with the zeros of the polynomial P n . We 
now have a closer look at the relation between the type I and type II multiple orthogonal 
polynomials and the eigenvectors of L n . The first n relations of the recurrence (2.2) can 
be written as 



( P (x) \ ( \ ( P (x 



\Pn-l(x) J 



+ 







\ 



Pn-2(x] 
\Pn-l(x) J 



(2.9) 



Substituting x = x^ n , £ = I, ... ,n, we then conclude that (Po(xe, n ) • • • Pn-i(xe, n )) T is the 
(only possible) right eigenvector of L n corresponding to the eigenvalue x^ n , normalized so 
that the first component is equal to 1. Note that the first component of a right eigenvector 
of L n can not be 0, otherwise each component would be 0. 

For the type I polynomials we can find something similar. Define for i = 1, . . . ,r the 



6 



vector polynomials 

2-1 

B n l \x) = xA n _ i+1 (x) - An-i(x) - J2 a n+j-i,j A n _ i+j+1 (x) (2.10) 

3=0 

r 

= a n-i+j,j (2-11) 
j=i 

where we use definition (2.11), with ag lt e 2 — 1, < £\ < £2 < r, in the cases 1 < n < i < r. 
The first n relations of the recurrence (2.4) then are 

(l^x) ■■■ A n (x)^j L n + (o_^_0 BP(x) ••• ^'(i)) 

n—r* 

= x (A^x) ■■■ A n (x)^j, (2.12) 

where r* = min(r, n). 

Definition 2.3 Construct the polynomials B$, i — 1, . . . , r, as in (2.10) and define for 
n eN and i = 1, . . . , min(r, n) the polynomials 

Q&W = det(sW(x) ••• ^(x) A,(x) S^)(x) ••• flW(x)), (2.13) 
fc G N. We also define 

B n (x) = det(s«(x) ■ ■ ■ B n r) (x)^j , neN. (2.14) 



Note that each linear combination of the components of the vector polynomials satisfies 
a similar relation as in (2.12). In particular, for the polynomials Q^ n , i — 1, . . . , min(r, n), 
we get 

(<?&(*) • • • Q<Sn(*)) L n + B n (x) e n T -i + l = * (Q { Sn(x) ■ ■ ■ Q$>n(*)) > ( 2 -15) 

where eg, 1 < £ < n is the £th unit vector in IR ra . In Lemma 2.4 below we will prove that 
the polynomial B n is equal to P n up to a multiplicative nonzero constant. Equation (2.15) 
then implies that the vectors 

(Qjfe) • • • Q?n(^)) T > * = !> • ■ • ,min(r,n), (2.16) 

are left eigenvectors of the matrix L n corresponding to the eigenvalue xg jn , if they are 
different from 0. 

Lemma 2.4 Let A n the vector polynomials defined by the recurrence relation (2.4) and 
its initial conditions (with Ajj 7^ ; j = 1, . . . ,r). Here we assume that the recurrence 
coefficients a^ T , £ > r, are different from 0. For the polynomials B n we then have that 
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B n (x) = 7„ P n (x), where the P n are the monic polynomials defined by the recurrence 
relation (2.2) and 



XV- A - ■ 

ln = iig*. h3 (-i)LfJ + LVJ, n = mr + S) 0<s<r. (2.17) 
Ue= r a e,r 



Proof. First of all we prove that B n is a polynomial of exactly degree n with leading 
coefficient j n . Using (2.11) we find that 



B n (x) = det(^B^(x) ••• B%\x) 

n+r-l , ^ ^ 

= J] a e>r detl A n+r {x) ■■■ A n+1 {x)\ (2.18) 

n+r— 1 

= II a *,r Yl si S n ( 7r ) A.+r,w(i)(a;) ■ ■ • A n+1;7T{r) (x), (2.19) 



where we denote by S r the set of permutations of r elements. Note that from the recurrence 
relation (2.4) and its initial conditions, with Ajj ^ 0, j — 1, . . . , r and a^ r ^ 0, £ > r, we 
obtain that 



A mr+8 , a (s) = As ' s x m + O (x™- 1 ) , l<s<r, m6NU{0}. (2.20) 

Combining this and (2.19) we see that, with n = mr + s, < s < r, 

deg(B n (x)) = deg(A n+r , s (x) . . . A n+r _ s+1>1 (x) A n+r _ S;r (x) . . . A n+1;S+1 (x)) = n, 



where we note that all the other terms in (2.19) have lower degree. The sign of the 
permutation corresponding to this term is (—1) Lf J + L^J so that j n is the leading coefficient 
of B n . 

To complete the proof we show that the polynomials satisfy the recurrence relation 
(2.2). Note that 

7n (-1)(-^ 



In- j U7=n-i a l 



n-i „ > l<j<", (2-21) 



r 
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where we assume ao, r = • • • = a r -\. r = 1. By (2.21) and (2.18), we then obtain that 

( B n+1 (x) ,B n (x) r B n-j(x) 
In + K,o - x) + 2^ a nJ J - 

V ln+1 In j = i In-j 

= (-ly 1 a n , r B n+l {x) + (a n , - x)B n (x) + £(-1)^" a nj -^-< f 



j=l Lle=n-j a t,r 

= (-l) r_1 Y[a e>r detl A n+r+1 (:r) ••• A n+2 (a;)) 

n+r-l , ^ ^ s 

+ (a„, - X) II °V det ( ^n+r(x) ■ ■ ■ A n+1 (x) J 

+ ]T(-l) (r - 1)j a nj J] a e , r det(l n _ j+r (x) ■ ■ ■ A n _ j+1 (x)). (2.22) 

3=1 l=n V J 

Form the recurrence relation (2.4) we first of all get that 
(-l) r_1 a n+r . :r det^4„ +r+1 (x) ■■■ A n+2 (x^j 

= (-l) r_1 det(-A n (x) - (a nfi - x)A n+1 (x) A n+r (x) ■ ■ ■ A n+2 (x)^j 

= -det(A n+r (x) ■■■ A n+2 (x) A n (x)^j - (a nfi - x) det (A n+r (x) ■■■ A n+1 (x)^j, 

so that the first two terms of (2.22) reduce to — 1 a-i,r det(A n+r (x) • • • A n+2 (x) A n {x)) . 

We apply this argument several times, knowing that, for j = 1, . . . , r — 1, 



a n+r , r dei(A n+r+1 - j (x) ■■■ A n+2 (x) A n (x) ■■■ A n _ j+1 (x) 
= (-l) r det^+^x) • • • A n+2 (x) A n (x) ■ ■ ■ A^^x^j 
+ (— l) r 3 a n j detl A n+r _j(x) ■ ■ ■ A n _j + i{x 



Finally we obtain that the expression in (2.22) is equal to 0. This proves the lemma. □ 
2.3 Generalized Christoffel-Darboux identity 



For orthogonal polynomials it is known that they satisfy the Christoffel-Darboux formula, 
see, e.g., [9]. In [16] the authors proved a generalized Christoffel-Darboux identity for ma- 
trix orthogonality of vector polynomials. This includes a generalized Christoffel-Darboux 
identity for the multiple orthogonal polynomials of type I and type II in the sense that one 
of the vector polynomials then has just one component. Here we show that this identity 
can be found as a natural consequence of Section 2.2. 
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Theorem 2.5 Suppose that the measures /ii, . . . ,/i r form a weakly complete system. For 
the corresponding multiple orthogonal polynomials we have 



(x - y) ]T P^ x (x)Q® n (y) = 7n(Pn(x)P n ^(y) - P n {y)Pn-i{x)) (2.23) 
k=i v 1 

and 

E ft-i(x)Qg,(x) = 7 n [P^P^ix) - P B (x)^_i(x)), (2.24) 
fc=i v 7 

i = l,... ,min(r, n), with Q$ n as in Definition 2.3 and 7„ as in (2.17). 

Remark 2.6 The Christoffel-Darboux formulas in Theorem 2.5 are related to the ge- 
neralized Christoffel-Darboux identity [16, (21)], with one vector polynomial having just 
1 component. In particular, for each i, expression (2.23) is a linear combination of the 
vector components in [16, (21)]. 

Proof. Define P{x) = (P (x) ••• P n _!(z)) T and $®(x) = {Q^ n {x) ■■■ Q$ n (x)) T , 
i — 1, . . . ,min(r, n). The equations (2.9) and (2.15) then reduce to 

L n P(x) + P n (x)e n = xP{x) (2.25) 

and 

Q^{x) T L n + B n (x)e n T - i+ i = xQ^(x) T . (2.26) 

Now multiply (2.25) on the right by the vector Q^\y) T and (2.26), evaluated at y, on the 
left by P(x), so that 

L n P(x)Q^(y) T + P n (x)e n Q^(y) T = xP(x)Q (i \y) T , 
P(x)Q^(y) T L n + B n (y)P(x)e n T - i+ i = yP(x)Q^(y) T . 

Next, take the trace of these two equations and subtract. Since for two arbitrary squared 
matrices C and D the property tr(CD) = tr(DC) holds, we then obtain 

(x - y) tr (P{x)Q^\y) T ) = P n {x)Q%{y) - B n (y)P n ^(x). (2.27) 

First of all we note that B n (y) = j n P n {y) by Lemma 2.4. Secondly, using (2.10) and 
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(2.11), we get 

= det(fl< 1 >(y) ••• Bt^iy) AM B% +1 \y) ••• flW(y)) 

= det( -A n _i(y) ••• -A n _ i+1 (?/) A n (?/) a„ +r _j_i ir A n+r _i(y) ••• a ni7 .A n+1 (y)) 

V v ' S v '/ 

2—1 r— 2 

rt+r— i— 1 , s 

= I] a e , r (-l)^Uet(A l _ i+r (y)..-A n _ i+1 (y)) 
= (-l) (r - 1)l f Tl o,, r ] B n ^y). 

\e=n-i J 

Combining this with Lemma 2.4 and (2.21) we can conclude that Q$ n (y) = 7n P n -i(y)- 
So, from (2.27) we finally obtain expression (2.23). If we let y — > x we also find (2.24). □ 



3 Multiple Gaussian quadrature 



Suppose hi, . . . , n r are measures for which all the moments exist. In this section we want 
to approximate a set of integrals of the kind 

/ fj(x) dfij(x), j = l,...,r. 

An obvious possibility is to approximate each of these integrals separately using Gaussian 
quadrature. However, when f 1 — --- — f r — fwe can consider to take the same set of 
quadrature nodes in each weighted quadrature formula. Our approximation then looks 
like 

f(x)d H (x) = Y,w < f,lf(xt,n) + Ei>\f), j = l,...,r, (3.1) 

for different quadrature nodes Xi >n , . . . , which is due to C. F. Borges [5]. So, we will 
not necessarily have the maximal order for each of the weighted quadrature formulas. 
However, the advantage of this choice is that we only need n evaluations of the function / 
instead of rn. In accordance with the case r = 1 we say that the set of weighted quadrature 
formulas (3.1) has vector order d = (di, . . . , d r ) E Nq if E^\f) = for each polynomial / 
of degree less than or equal to dj, j — 1, . . . , r. 

Our goal is to maximize the vector order of (3.1) along the proper multi-indices. A nec- 
essary requirement is then that each of the quadrature formulas is interpolating, which 
means that we have vector order at least (n — l,...,n — 1). This corresponds to the 
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conditions 



W S=/ h,n(x) dfij(x), £=l,...,n, j = l,...,r, (3.2) 
where we denote by 

h,n( x )= II ZT^~' £=l,...,n, 

the fundamental polynomials of Lagrange interpolation. The only freedom we still have 
then consists of the choice of the set of quadrature nodes. In [5] C. F. Borges already 
proved that the maximal vector order (along the proper multi-indices) is obtained if we 
choose the quadrature nodes to be the zeros of the type II multiple orthogonal polynomial 
P n , corresponding to the measures /ii, . . . , \± r . We recall this theorem and give a proof for 
completeness. Of course we need that the proper multi-index V n is normal so that P n is 
uniquely defined and has exact degree. Furthermore, we need that the zeros of P n are 
simple. Note that, in the case of positive measures, these conditions are satisfied if the 
measures form, for example, an AT system or Angelesco system [15,17]. 

Theorem 3.1 (Borges) Suppose that the measures //i,...,/x r form a weakly complete 
system and that P n , the type II multiple orthogonal polynomial of degree n, has simple 
zeros. We then speak of multiple Gaussian quadrature if the weighted quadrature formulas 
in (3.1) are interpolating and the quadrature nodes are the zeros of P n . In this case we 
have vector order (n — l)e + V n , where e— (1, . . . , 1) G N r . Furthermore, these quadrature 
formulas do not have vector order (n — l)e + v n +i- 

Proof. We prove that, with the conditions of the theorem, the j'th weighted quadrature 
formula in (3.1) has order n — 1 + v n {j)-> 3 — 1, • • • ,r. Let hj be a polynomial of degree 
at least n and at most n — 1 + v n {j)- (If v n (j) = 0, there is nothing to prove.) Denote by 
T^li the interpolating polynomial of hj at the zeros xi >n , . . . , x n>n of the type II multiple 
orthogonal polynomial P n . The polynomial hj — can then be written as 

h J (x)-T^ 1 (x) = P n (x)R^(x), 

where is a polynomial of degree at most u n (j) — l. Since the jth weighted quadrature 
formula is interpolating, we have 

/n ,. n 

hj(x) dHj{x) - ]T hj(xe tn ) = / hj( X ) dfljix) - J^n-lfon) W g 

e=i J t=\ 

= ( (hj^-T^xj) dvj{x) 



j P n {x)RP\x) dfij{x). 



By the orthogonality conditions of the type II polynomial P n we then see that this is 
equal to 0. 
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The x^ n are the zeros of the polynomial P n , so 

n r 
J2 P n(xe,n) A n+l,j( X e,n) «fi = 0. 

e=i 3=1 

Since P n +i is a normal index, this is different from f P n (x)Y?j=iA n +ij(x) d/j,j(x), which 
means that we do not have vector order (n — l)e + z4+i- □ 

We now study the theory of multiple Gaussian quadrature from the practical point of view. 
In particular we give a link with the eigenvalue problem of the banded lower Hessenberg 
matrix L n , see (2.8). First of all, from Section 2.2 we know that the zeros of the polynomial 
P n are the eigenvalues of L n . In the theorem below we show that, for multiple Gaussian 
quadrature, the quadrature weights can be expressed in terms of the corresponding left and 
right eigenvectors. This extends the expressions found in the case of Gaussian quadrature, 
see e.g. [10, Chapter 3, §2.3] and [11]. 

Theorem 3.2 Assume that the measures /xi, . . . , fj, r form a weakly complete system and 
that the type II multiple orthogonal polynomial P n has simple zeros x it7l , . . . , x n>n . Let v^ n 
be the right eigenvector of the matrix L n corresponding to the eigenvalue x^ n , with first 
component equal to 1. Similarly, denote by ue t n the corresponding left eigenvector with the 
keth component, the first non-zero component, equal to 1. Here 1 < h < min(r, n). In the 
case of multiple Gaussian quadrature we then have, for £ = 1, . . . , n, 

-y I m i n (ii n ) \ 

w^ n = _>rp - J2 C hk u e , n (k) , j = l,...,r, (3.3) 

U £,n V e,n \ k=1 J 

where 

C jtk = I P k -i(x) d(ij(x) (3.4) 
= V(-l) fe+i mg?! 7 * , 1< k < j < r. 

Here the constants Cjj, j = 1, . . . , r are different from and wf^ — ; j — 1, . . . , ki — 1. 
Proof. In this proof we fix £ E {1, ... , n}. From Section 2.2 we know that 

By Remark 2.1 we have a m + r ~i,r 7^ 0, m > 1. Then it is clear that at least one of the 
first min(r, n) components of the left eigenvector ui^ n and at least one of the last min(r, n) 
components of v^ n is different from 0. So, there exists an i£ e {1, . . . , min(r, n)} for which 
P n -i e ( x i,n) 7^ 0. Taking x = X£ >n in (2.24) with i = i£, we then find 



J2 P k-l(xe,n)Qk!l( X i,n) = lnP' n { X i,n)Pn-H{x£,n) + 0, (3.5) 
k=l 
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because x^ n is a simple zero of P n . Consequently, the vector (2.16) with i = in is different 
from and 



UP; 



... Q 1 Q£|l,n(^,n) ... Qgfo,n) \ T 



fef-i 

Applying (3.5) we then find that 

7nPn( X e,n)Pn-i e (xe,n) 



(3.6) 



(3.7) 



Next, since the weighted quadrature formulas are assumed to be interpolating, we obtain 
from the Christoffel-Darboux formula (2.23) with % — %(_ and y = x^ n that 



™el = I kn(x) dfij(x) 



Pn(x) 



d/ij(x) 



1 n f 

Combining this with (3.7) and using the orthogonality conditions of the type II polynomial 
P n we then finally get 



7/ ,0') _ _ 

w l,n ~ r+T 



u, 



E 



in) 



P k -i(x) dnj(x), j = l,...,r, 



which proves (3.3). Finally we note that the proper multi-indices are normal, so we have 
that Cjj 7^ 0, j — 1, . . . , r. 



□ 



Remark 3.3 There exists a one to one mapping between the sets of constants Cj^ G C, 
1 < k < j < r, (with Cjj ^ 0, 1 < j < r) and A i:j e C, 1 < j < i < r, (with A jd ^ 0, 
1 < j < r). By the orthogonality conditions of the multiple orthogonal polynomials of 
type I and type II we obtain that 



j=k 



J2 A i,j °j,k = / Pk-l(x) J2 A iJ d Vj( X ) 



3=1 



0, k — 1, . . . , i — 1, 

1 , k = i, 



i = l,...,r. 



In particular, each of these sets can be found from the other one by solving the linear 

systems 



^1,1 ^2,1 ' ' ' <^i,l 
6*2,2 ' ' " Ci,2 

o ••• o a* 



Ai,2 





V 1 / 



i = l,...,r, 



(3.8) 
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and 



Ai,i • 



V- 



( a, \ 



a 



\ C rt i ) \l)J 



i = l,...,r, 



(3.9) 



respectively. 



4 Conclusion 



In Theorem 3.2 we proved that in the case of multiple Gaussian quadrature the quadrature 
nodes and weights in (3.1) can be found explicitly by solving the left and right eigenvalue 
problem of a banded lower Hessenberg matrix L n . Here L n contains the recurrence co- 
efficients of the multiple orthogonal polynomials corresponding to the set of measures 
/ii, . . . , fi r as in (2.8). Further, we only need a set of constants Cj^ G C, 1 < k < j < r, 
(with Cjj 7^ 0, 1 < j < r) which can be found from the initial values of the type I multiple 
orthogonal polynomials by solving the linear systems (3.8). 

By Theorem 3.1 the weighted quadrature formulas preserve the orthogonality conditions 
(and normalization) for the finite set of polynomials A 1 , . . . , A n and P , Pi, ... , P n -i- So, 
these are also the type I and type II multiple orthogonal polynomials corresponding to 
the set of discrete measures 

n 

^3,n = X! W £,n^x e , n , j = 1 , ■ ■ ■ , r , 

£=1 

with finite support. Then note that Remark 3.3 is a nice illustration of Remark 2.2. Since 
we explained how to find these discrete measures starting from the recurrence coefficients 
and the set of initial values A it j e C, 1 < j < % < r, (with Ajj ^ 0, 1 < j < r), this gives 
rise to a constructive proof of the spectral theorem for a finite set of multiple orthogonal 
polynomials. In the case of an infinite set some results for the spectral theorem were 
already obtained in [13,16]. However, finding necessary and sufficient conditions on the 
recurrence coefficients to have positive orthogonality measures on the real axis is still an 
open problem. 
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